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ABSTRACT 


Integrated radio-spectrum of Cas A in continuum was analyzed with special 
emphasis on possible high frequency spectral curvature. We conclude that the 
most probable scenario is that Planck’s new data reveal the imprint of non-linear 
particle acceleration in the case of this young Galactic supernova remnant (SNR). 

Subject headings: ISM: individual (Cas A) - radiation mechanisms: non-thermal - 
acceleration of particles - ISM: supernova remnants - radio continuum: ISM 
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Introduction 


Planc!^ is the third generation space mission to measure the anisotropy of the cosmic 
microwave background. It observed the sky in nine frequency bands covering 30 - 857 
GHz with high sensitivity and angular resolution from 5' to 31' (Arnaud et ah 2014 and 
references therein). The Low Frequency Instrument covers the 30, 44, and 70 GHz bands 
while the High Frequency Instrument covers the 100, 143, 217, 353, 545, and 857 GHz 
bands. Planck’s sensitivity, angular resolution, and frequency coverage make it a powerful 
instrument for Galactic and extragalactic astrophysics as well as cosmology. 


Gassiopeia A (Gas A) is a very bright and young supernova remnant (SNR), likely due 
to a historical supernova around 353 — 343 yr ago (Ashworth 1980; Fesen et al. 2006). Its 
estimated distance is 3.33 ± 0.10 kpc (Alarie et al. 2014). Recently, this SNR was observed 
by Planck and the results were published in Arnaud et al. (2014). It was shown that Gas 
A is a distinct compact source from 30 - 353 GHz but becomes confused with unrelated 
Galactic clouds at the highest P/anc/c frequencies (545 and 857 GHz). The apparent excess 
radiation at 217 and 353 GHz is proposed to be due to coincidental peak in the unrelated 
foreground emission or to cool dust in the supernova remnant (Arnaud et al. 2014). 


The slightly concave-up forms of radio-spectra were detected for some young SNRs 
(Reynolds & Ellison 1992). The main reason for high frequency concave-up curvature 
in radio-spectra of young SNRs should originate in non-linear diffuse shock particle 
acceleration (see Urosevic 2014 and references therein). Due to positive identihcation of 


^Planck is a project of the European Space Agency (ESA) with instruments provided 
by two scientific consortia funded by ESA member states, with contributions from NASA 
(USA) and telescope reflectors provided by a collaboration between ESA and a scientific 
consortium led and funded by Denmark. 
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infrared synchrotron radiation from Cas A, Jones et al. (2003) gave indication that the 
radio-spectrum of this SNR should be concave-up. 

This paper focuses on the analysis of origin of the high-frequency curvature in 
radio-spectrum of Cas A. 


2. Analysis and Results 

There are various observations of SNR Cas A (Green 2014). Flux densities at different 
frequencies for SNR Cas A were taken from different papers. Different samples for data 
analysis were used to account for differences in data gathered from various literature. The 
hrst sample includes data from Arnaud et al. (2014), Baars et al. (1977) and references 
therein as well as from Mason et al. (1999). The original data (data for original epochs) 
are taken from Table 2 of Baars et al. (1977), not the scaled ones. The second sample is 
formed by addition of data from Hurley-Walker et al. (2009) and Liszt & Lucas (1999). The 
frequency range for both samples is 0.55 - 353 GHz. The integrated flux density of 52 ± 7 
Jy at 353 GHz was adopted in accordance with the discussion presented in Arnaud et al. 
(2014). 

Before actual analysis, the flux densities were appropriately scaled to account for 
secular fading. The epoch (2010.0) is chosen to coincide with Planck’s intermediate 
astrophysics results (Arnaud et al. 2014). Two scaling relations were used in this paper. 
The first one, commonly used, is taken from Baars et al. (1977) in the following form 

d{u) = a- 61ogi/[GHz], 

Ad{i/) = Aa + Ab \ log i^[ghz]I, 

a = (0.0097 ± 0.0004), b = (0.003 ± 0.0004), (1) 

where d{v) is the secular decrease in the flux density at given frequency and Ad{v) is an 
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appropriate error estimate. 

On the other hand, it was debated that this scaling relation is not appropriate 
(O’Sullivan & Green 1999) especially for the lowest radio-frequencies (Helmboldt & Kassim 
2009 and references therein). 

Due to the lack of accuracy of the above mentioned relation, the scaling proposed in 
Vinyaikm (2014) was also used in the following form 

d{iy) = a-b\n z/[ghz] - c 
Ad{iy) = Aa -h A6 I In z/[ghz] I + Ac z^[ghI]> 
a = (0.0063 ± 0.0002), b = (0.0004 ± 0.0001), 

c = (0.0151 ± 0.0016) X 10"^ (2) 


The results obtained analyzing data formed by these two scaling relation are compared. 
The flux density for the desired epoch and the appropriate error estimate are calculated 


St^{iy) = St^iiy) (1 - d{u)) , T = t2 - h, 

ASM = SM 


ASt,{iy) Ad{u) 


SM 


1 — d{v) 


( 3 ) 


It is worth noting that radio-spectrum of Gas A shows low-frequency cut-off due to 
thermal absorption (Vinyaikm 2014) or possibly synchrotron self-absorption. As the main 
interest of this paper involves analysis of high frequency part of the Gas A radio-spectrum, 
hrst, only data at frequencies higher than around 550 MHz were analyzed (for both scaling 
relation). 

To show that the observed high-frequency radio-spectrum of Gas A is actually curved, 
the principal component analysis (PGA) was first applied (Babu & Feigelson 1996). The 
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Fig. 1.— The results of PCA on data from the first sample. Scaling relation taken from 
Baars et al. (1977) was used. PCI and PC2 correspond to the hrst and second principal 
component, respectively. In the upper left figure are the results of PCA on the radio-data 
between 550 MHz and 30 GHz. Departure from the linear relationship is obvious when the 
observation at 30 GHz is added (upper right hgure), both data at 30 and 32 GHz are added 
(down left hgure) and when the whole range of frequencies up to 353 GHz is used (down 
right hgure) 
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principal components are the eigenvectors of the covariance matrix. Original data are 
represented in the basis formed by these vectors. The hrst principal component accounts 
for as much of the variability in the data as possible, i.e. it has the largest possible variance. 
The second principal component has the highest variance possible under the constraint 
that it is orthogonal to the hrst one. In that sense, PCA includes calculation of the highest 
data variability direction which can be used to detect departures from the pure power-law 
spectra (linear £t in log-log scale). 

In Figure 1, the results of PCA on logarithmically transformed data (logi^, logFi,) from 
the hrst sample are presented. PCI and PC2 correspond to the hrst and second principal 
component, respectively. The plotted values represent the original data in the basis of 
principal components (zero centered data multiplied by the rotation matrix whose columns 
contain the eigenvectors). The covariance of these values represents a diagonal matrix with 
squares of standard deviations of the principal components eigenvalues as its elements. In 
the upper left hgure are the results of PCA on the radio-data between 550 MHz and 30 
GHz. It is clear that the radio-spectrum follows pure power-law relation. Addition of the 
observation at 30 GHz (upper right hgure), both data at 30 and 32 GHz (down left hgure) 
and if whole range of frequencies up to 353 GHz is analyzed (down right hgure), shows that 
the radio-spectrum of Gas A is positively curved (concave up) above around 30 GHz. The 
similar results are obtained for the radio-data scaled by the relation of Vinyaikm (2014). 

If data taken from the second sample are used, similar results are obtained although the 
scatter is more prominent. The same conclusions hold if both scaling relations are used. 

As the high-frequency radio-spectrum of Cas A is clearly curved, it is not appropriate 
to be fitted by the pure power-law. In the most simple case, it can be fitted by the 
two power-laws, below and above 30 GHz (excluding, of course, the low-frequency data 
inhuenced by the low-frequency cut-oh). On the other hand, it is known that one of the 
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Table 1: Varying power-law best fitting parameters (0.55 - 353 GHz). 


sample (scaling) 


a 



a 


( 

[k) 

R2 

1 (Baars at al. 1977) 

0.790 

± 

0.016 

0.056 

± 

0.008 

23.11 

(34) 

0.997 

2 (Baars at al. 1977) 

0.806 

± 

0.016 

0.062 

± 

0.008 

38.55 

(42) 

0.995 

1 (Vinyaikm 2014) 

0.841 

± 

0.014 

0.065 

± 

0.007 

27.33 

(34) 

0.997 

2 (Vinyaikm 2014) 

0.854 

± 

0.013 

0.069 

± 

0.007 

50.54 

(42) 

0.995 


repercussions of non-linear effects of particle acceleration is the curved particle spectrum 
which in turn gives rise to the curved synchrotron spectrum (Ellison & Eichler 1984). 

In that sense, it is more appropriate to represent radio-spectrum of young SNRs with 
approximative relation that represents varying power-law of the form 

%!(>') = %l(lGHz) ('I) 


where a is the standard synchrotron spectral index and a is the parameter of spectral 
curvature which should be positive due to the non-linear behavior of DSA (Houck & Allen 
2006; Vinyaikm 2014 and references therein). Of course, more general relation applied to 
the larger range of energies would include appropriate high-frequency synchrotron spectral 
cut-off which we do not take into account in this analysis. 


The results of weighted least-square hts to data for this model applied to both samples 
and for different scaling relations are summarized in Table 1. Fit properties are given in 
form of and dehned as 

N 


= ^Wi {vi - fixi] (3i ,..., /3p)f 


2=1 


1 

) 2 ’ 


= 1 - 


X 


TSS’ 


TSS = '^Wi{yi-yf , y= 

Li=i uy 


2=1 


( 5 ) 
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In above relations, / is the predicted value from the £t, f5j are fit parameters {j = 1, ...,p) 
and y is the weighted mean of the observed data y* at particular Xi, i = 1,N. Wi are 
the weights applied to each data point. TSS is the so called total sum of squares. It is 
worth stressing that in the case of non-linear models, the number of degrees of freedom is 
generally unknown, i.e. it is not possible to compute the value of reduced or adjusted 
(Andrae et al. 2010). In Table 1, parameter k = N — p is given for convenience. In 
non-linear models, k does not always represent the exact number of degrees of freedom 
(Andrae et al. 2010). 

There are no significant differences in the fit results for both samples and for different 
scaling relations. Again, the exclusion of the data point at 353 GHz does not change the 
results signihcantly. It must be noted that the smoothly curved spectral shape represented 
by equation (4) gives a better fit than the model that assumes two power laws. 

In Figure 2, the weighted least-square fit to the data from the hrst sample, for the 
varying spectral index model is presented. Scaling relation taken from Vinyaikm (2014) was 
used. Diamond symbols indicate Planck’s data. 


It is also appropriate to check if obtained values of above mentioned parameters, a and 
a, are consistent with the ones obtained when low-frequency cut-off is taken into account. 
To this end, the previous model was adjusted in the following way 


•S[Jy|(l') 


S|j,,(lGHz) 


-a+a log !/[ghz] 
^[GHz] 


^-rou 


- 2.1 

[GHz] ^ 


( 6 ) 


and is of the same form as the one stated in Vinyaikm (2014). tq represents the optical 
depth at 1 GHz. 

For this analysis new data samples were formed and analyzed. The third sample 
includes data from Tables 3 and 4 from Vinyaikm (2014) for the epoch 2015.5 as well as 
the Planck’s data from Arnaud et al. (2014) appropriately scaled to match the same epoch. 
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The frequency range for this sample is 0.0056 - 353 GHz. The fourth sample includes all 
the data from Table 2 presented in Baars et ah (1977) with addition of data from Arnaud 
et al. (2014), Helmboldt & Kassim (2009) and Mason et al. (1999). Finally, the hfth sample 
incorporates the fourth one with the addition of data from Hurley-Walker et al. (2009) and 
Liszt & Lucas (1999). Data in fourth and hfth sample are scaled to the epoch of Planck’s 
observations and the frequency range for these samples is 0.01005 - 353 GHz. In all these 
cases, scaling relation from Vinyaikm (2014) was only used. 

The results of weighted least-square hts to data in this case are summarized in Table 
2. There are no signihcant differences in the ht results for all samples. In Figure 3 the 
weighted least-square ht to the data from the fourth sample, for the varying spectral index 
model with low-frequency cut-oh is presented. Diamond symbols indicate Planck’s data. 

The results are not signihcantly diherent than estimated in Vinyaikm (2014). On 
the other hand, the estimated values for spectral index around 1 GHz as well as for the 
curvature parameter are less than the ones obtained using the relation (4) on high frequency 
part of the radio-spectra alone. 

The apparent excess radiation at 353 GHz, discussed in Arnaud et al. (2014), is also 
present in this analysis (see Figure 3). It may be due to the presence of cool dust in 
northern and western parts of Gas A (Dunne et al. 2009; Vinyaikm 2014). It should be 
noted that images at 600 GHz observed with the Herschel Space Observatory (Barlow 

Table 2. Varying power-law with low-frequency cut-off best htting parameters. 


sample 

frequency range [GHz] 

a 

a 

"ro 

(fc) 


3 

0.0056 - 353 

0.762 ± 0.005 

0.033 ± 0.003 

(7.738 ± 0.574) x 10“® 

83.05 (29) 

0.994 

4 

0.01005 - 353 

0.760 ± 0.005 

0.025 ± 0.003 

(8.968 ± 0.869) x 10“® 

75.79 (46) 

0.994 

5 

0.01005 - 353 

0.765 ± 0.005 

0.026 ± 0.003 

(9.209 ± 0.875) x 10“® 

111.03 (54) 

0.991 



10 ^ 


10 ' I_ 

0.001 


a = 0.760 + 0.005 
a = 0.025 + 0.003 
To = (8.968 + 0.869)x10"=’ 


0.010 0.100 1.000 

V [GHz] 



10.000 100.000 1000.000 


Fig. 3.— The weighted least-square ht to the data from the fourth sample, for the varying 
spectral index model with low-frequency cut-off. Diamond symbols indicate Planck’s data 
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et al. 2010) at much higher angular resolution than Planck (6" - 37") show that the 
non-synchrotron microwave emission may be a combination of both cold interstellar dust 
and freshly formed dust. Arnaud et al. (2014) stressed that this excess could potentially 
be due to a coincidental peak in the unrelated foreground emission or to cool dust in the 
supernova remnant, which is marginally resolved by Planck. 

The spinning-dust emission (electric dipole radiation from the rapidly rotating dust 
grains) could generally contribute signihcantly at high radio-frequencies, especially around 
10 - 100 GHz (Draine & Lazarian 1998ab; Ali-Hai'moud et al. 2009; Scaife et al. 2007; 
Silsbee et al. 2011; Stevenson 2014). On the other hand, in the case of Gas A this emission 
component is negligible. 

Finally, Vinyaikin (2014) showed that the observed slowing of the secular variations of 
the radio flux density of Gas A with decreasing frequency at decimeter wavelengths can be 
explained by a decrease in the optical depth of a remnant H ll zone around Gas A with time 
due to recombination of hydrogen atoms. Due to the presence of the H ll region it may be 
interesting to fit the radio spectrum of Gas A with the following expression 


5(1.) = 5i z/-“ e-"-' + 52 (1 - e-"") , 


( 7 ) 


This formula incorporates both synchrotron radiation as well as the thermal absorption 
and thermal bremsstrahlung emission from the zone in front of the synchrotron emitting 
region. Of course, this is a rather naive model. Supernova remnants are 3D structures, the 
fact that this relation does not take into account. On the other hand, it can serve as an 
approximative model for qualitative analysis. 

If frequencies are in GHz then the (non-thermal) synchrotron flux density at 1 GHz is 
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simply 5'^"'"(lGHz) = Si but for the thermal component holds 

5^(lGHz) = ^2 (1 - e-"°) , 

AS’'(lGHz) = SyiGHz) + At„ j . (8) 

SO that the contribution of thermal emission in integral radiation at 1 GHz is given by 

^ _ S^ilGRz) 

^ ~ ^t(igHz) + ^nt(igHz)’ 

( AST(igHz) I A5NT(igHz) 

ST(igHz) + SNT(igHz) 

^ , SNT(igHz) 

5T(1GHz) 

The results of weighted least-square hts to data in this case are summarized in Table 3. 
Again, there are no signihcant differences in the £t results for all samples. The weighted 
least-square fit to the data from the fourth sample, for the test-particle synchrotron model 
with thermal absorption and emission is presented in Figure 4. Diamond symbols indicate 
Planck’s data. 

The presence of thermal emission from the Hll region associated with this SNR is 
estimated to be around 1-2 % in the total flux density at 1 GHz. Addition of this component 
significantly improves fit to the data but discards the curvature due to non-linear effects of 
particle acceleration. It is worth noting that addition of the curvature parameter a to the 
previous model, i.e. varying power-law model with thermal absorption and emission in front 

Table 3. Test-particle synchrotron model with thermal absorption and emission best 

fitting parameters. 



sample 

frequency range [GHz] 

OL 

"^0 

? [%] 

(fc) 

B? 

3 

0.0056 - 353 

0.769 ± 0.006 

(5.211 ± 0.341) X 10“® 

1.653 ±0.005 

53.97 (29) 

0.996 

4 

0.01005 - 353 

0.785 ± 0.008 

(7.846 ± 0.745) x 10“® 

1.913 ±0.008 

47.35 (46) 

0.996 

5 

0.01005 - 353 

0.788 ± 0.008 

(7.535 ± 0.624) x 10“® 

1.887 ± 0.007 

74.12 (54) 

0.994 




Fig. 4.— The weighted least-square ht to the data from the fourth sample, for the test- 
particle synchrotron model with thermal absorption and emission. Diamond symbols indicate 
Planck’s data 




16 


of synchrotron emitting region, does not produce acceptable fits to all samples. Negative, 
i.e. non-physical in sense of the non-linear DSA effect^, values for the curvature parameter 
a are preferred. Bounding the parameters to appropriate physical ranges of values do not 
improve hts that are also highly sensitive to initial guesses for the parameters. On the 
other hand, simple model represented by equation (7) hts to data very well (see Figure 4). 
The hux density estimate at 353 GHz is within the error bars of the measured value at this 
frequency for data samples 4 and 5 (see Figure 4). 


3. Discussion 

It is known that some young SNRs, such as Gas A and Kepler exhibit evidence of radio 
spectral variations from one region to another (Anderson & Rudnick 1996; DeLaney et ah 
2002). This can lead to a curved composite spectrum even if the synchrotron spectrum for 
each region is just a power-law (Allen et ah 2008). On the other hand, the results of Jones 
et ah (2003) indicate that the spectra of particular small features in Gas A hatten with 
increasing energy and Allen et al. (2008) found that their results can be explained with 
model that includes curvature of particle spectra of b = 0.06 ± 0.01. Gurvature of particle 
spectra b is in fact four times spectral curvature a so that our results, although slightly 

^When the pressure of cosmic-ray (GR) particles, produced at the shock wave, is included 
in the analysis, such a modihcation of shock structure implies changes in the spectrum of 
accelerated electrons. The energy spectrum of low energy electrons becomes softer (the 
radio spectrum is steeper) and the spectrum of high energy electrons becomes harder (the 
radio spectrum is shallower) than that of linear DSA. Due to signihcant GR production, 
the concave-up (positively curved) radio spectra is expected in the case of young SNRs (see 
Urosevic 2014, and references therein for more details). 
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higher, are not in disagreement with conclusions of Allen et ah (2008). 

For additional support to our claims, we used the same model for synchrotron radiation 
as in Allen et ah (2008), with one simplification that the non-thermal particle distribution 
function lacks high energy cut-off 

N(E)dE = K ( — j , (10) 

where b = 4a is the curvature of particle spectra, Eq = 1 GeV and F is particle spectral 
index at Eq. Flux density is a combination of synchrotron radiation from such an ensemble 
and thermal absorption at low frequencies in a similar manner as in equation (6). The 
results, given in Table 4, show that the simplified model based on the one from Allen et ah 
(2008) is consistent with our previous results. 

Finally, Dunne et ah (2009) showed that polarized sub-millimeter emission is associated 
with the SNR and that the excess polarized sub-millimeter flux at 353 GHz is due to cold 
dust within the remnant. They noted that there is no currently known way to produce such 
a polarized sub-millimeter emission from a synchrotron process. 

Keeping in mind the fact that inclusion of thermal bremsstrahlung emission (equation 
7) apparently leads to a better £t to data (Table 3), it must be noted that such an 
interpretation leads to a less contribution of dust emission at 353 GHz, which is in contrast 

Table 4. Simplified model based on the one from Allen et ah (2008) best fitting 

parameters. 


sample 

frequency range [GHz] 

a. 

a 

ro 

(fc) 


3 

0.0056 - 353 

0.775 ± 0.007 

0.032 ± 0.001 

(8.032 ± 0.500) X 10-® 

79.57 (28) 

0.994 

4 

0.01005 - 353 

0.739 ± 0.005 

0.025 ± 0.003 

(8.998 ± 0.867) x 10“® 

75.72 (45) 

0.994 

5 

0.01005 - 353 

0.777 ±0.008 

0.026 ± 0.003 

(9.268 ± 0.872) x 10“® 

110.90 (53) 

0.991 
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to the conclusion of Dunne et al. (2009). In that sense, we conclude that the more probable 
scenario is that Planck’s new data actually reveal the imprint of non-linear particle 
acceleration in the case of this young Galactic SNR. Of course, we can not fully dismiss the 
possible contamination by adjacent Hll region. 

The additional attempt to support the non-linear particle acceleration hypothesis may 
come from the fit using non-linear synchrotron model with thermal absorption and dust 
Planck-like emission (Figure 5 and Table 5). To that end, the Herschel’s total flux densities 
from Barlow et al. (2010) were added to data samples 3-5 forming samples 6-8. Again, 
similar results are obtained. Of course, this is a simplihed model which does not distinguish 
between different dust emission components. On the other hand, the spectral curvature due 
to the non-linear particle acceleration is clearly present in this case too (see Table 5). These 
values for a are also more in accordance with the result of Allen et al. (2008). 

It is also worth mentioning that Pohl et al. (2015) showed that stochastic re-acceleration 
of electrons downstream of the forward shock can explain the soft spectra observed from 
many Galactic SNRs. They also noted that, generally, interiors of the SNRs produce 
slightly steeper radio spectra than does the shell where re-acceleration occurs. 

Table 5. The non-linear synchrotron model with thermal absorption and dust emission 

best htting parameters. 


sample 

frequency range [GHz] 

a. 

a 

ro 

(fc) 


6 

0.0056-4285.714 

0.756 ±0.005 

0.022 ± 0.004 

(6.807 ± 0.579) x 10"® 

60.86 (33) 

0.995 

7 

0.01005 -4285.714 

0.760 ± 0.005 

0.020 ± 0.004 

(8.559 ± 0.860) x 10“® 

69.55 (50) 

0.994 

8 

0.01005 -4285.714 

0.764 ±0.005 

0.021 ±0.003 

(8.800 ± 0.865) X 10“® 

102.01 (58) 

0.992 




Fig. 5.— The weighted least-square ht to the data from the seventh sample, for the non¬ 
linear synchrotron model with thermal absorption and dust emission. Diamond symbols 
indicate Planck’s and triangles Herschel’s data 
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4. Conclusions 

The main conclusions that can be drown from this analysis are 

(1) Planck’s data support the observation that the radio continuum of SNR Cas A 
flattens with increasing frequency. The spectrum becomes positively curved above around 
30 GHz. 

(2) Alternative explanations of curvature are investigated: signihcant non-linear effects 
of particle acceleration, presence of thermal emission due to the emission of associated H ll 
region, and the presence of dust emission. 

(3) Different samples of data were analyzed to account for differences in data gathered 
from various literature. The presented results are not signihcantly sensitive to differences 
between samples. 

(4) The results presented in this paper agree with the conclusions stated in Vinyaikm 
(2014), i.e. the observed flattening of the radio spectrum at millimeter wavelengths is most 
likely primarily due to a flattening of the radio synchrotron spectrum of Cas A itself. In 
other words, non-linear effects of particle acceleration are possibly mainly responsible for 
the apparent high-frequency curvature in Cas A radio spectrum. 

We wish to thank the anonymous referee for useful suggestions which substantially 
improved this paper. This work is part of Project No. 176005 ’’Emission nebulae: structure 
and evolution” supported by the Ministry of Education, Science, and Technological 
Development of the Republic of Serbia. 
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